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ABSTRACT 

The idea that wc Uvc in a Universe undergoing a period of acceleration is a new. 
yet strongly held, notion in cosmology. As this can, potentially, be explained with 
a modification to General Relativity we look at current cosmological data with the 
purpose of testing aspects of gravity. Firstly we constrain a phenomenological model 
(mDGP) motivated by a possible extra dimension. This is characterised by a parameter 
a which interpolates between a = (LCDM) and a = 1 (the Dvali Gabadadze Porrati 
(DGP) model). In addition, we analyse general signatures of modified gravity given 
by the growth parameter 7 and power spectrum parameter E. We utilise large angular 
scale {9 > 30 arcminutes) Weak Lensing data (CFHTLS-wide) in order to work in 
the more linear regime and then add, in combination, Baryon Acoustic Oscillations 
(BAOs) and Supernovae. We subsequently show that current weak lensing data is not 
yet capable of constraining either model in isolation. However we demonstrate that 
even at present this probe is highly beneficial, for in combination with BAOs and 
Supernovae we obtain a < 0.58 and a < 0.91 at la and 2(7, respectively. Without 
the lensing data no constraint is possible. Inclusion of all angular scales {1^9^ 230 
arcminutes) allows a noticeable but not significant increase in the constraining power 
of the joint probes with a < 0.56 and a < 0.86 also at la and 2a, respectively. Both 
analyses correspond to a disfavouring of the fiat DGP braneworld model (a = 1) at 
over 2a. We highlight these are insensitive to potential systcmatics in the lensing data 
such as an underestimation of the shear at high rcdshift. For the growth signature 7 we 
show that, in combination, these probes do not yet have sufficient constraining power. 
Finally, we look beyond these present capabilities and demonstrate that Euclid, a 
future weak lensing survey, will deeply probe the nature of gravity. A Icr error of 0.104 
is found for a (/max = 500) whereas for the general modified signatures we forecast la 
errors of 0.045 for 7 and 0.25 for Sg (/max = 500), which is further tightened to 0.038 
for 7 and 0.069 for Sq (/max = 10000). 
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* Based on observations obtained with MegaPrime/MegaCam, 
a joint project of CFHT and CEA/DAPNIA, at the Canada- 
France-Hawaii Telescope (CFHT) which is operated by the Na- 
tional Research Council (NRC) of Canada, the Institut National 
des Sciences de I'Univers of the Centre National de la Recherche 
Scientifique (CNRS) of France, and the University of Hawaii. This 
work is based in part on data products produced at Tbraplx and 
the Canadian Astronomy Data Centre as part of the Canada- 
France-Hawaii Telescope Legacy Survey, a collaborative project 
of NRC and CNRS. 



1 INTRODUCTION 

General relativity, a cornerstone of physics, is arguably one 
of our greatest intellectual achievements. It is not only el- 
egant and physically motivated, but it makes a whole host 
of predictions including gravitational waves, the anomalous 
precession of Mercury and the deflection of light-all of which 
have been verified. 

Further still, due to its weak yet accumulative nature, 
gravitation is a force of the large scale. And so cosmology, 
being the study of the Universe, is a field performed within 
the formalism of this theory. It is in this way that Cosmology 
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itself becomes yet a further test for gravity. Considering this, 
today's cosmologists have been posed the most tantahsing 
problem: given that recent precision data from Supernovae, 
the Cosmi c Microwave Backg ro und and large s cale struc- 
ture fE.g. lAstier et al.1 (|200g ). iDunklev et (|2008l ) and 
IPercival et all lj2007l 'lV all indicate that the Universe is un- 
dergoing a period of cosmic acceleration, do we treat this as 
evidence that Einstein's theory of gravitation, or the manner 
in which we implement it, is incomplete? Or, more conven- 
tionally, should we stand by this successful theory and invoke 
some new matter component-Dark Energy-to explain it? 

Here we study the former and investigate the idea that 
General Relativity is not general enough. We do not at- 
tempt to motivate a new theory of gravity itself but instead 
aim towards testing existing theories and aspects of gen- 
eral theories with current and future data. In Section [2] we 
review the concept of m odified gravity, includ ing a gener- 
alised braneworld model (|Dvali fc Turneill2003l ) we go on to 
constrain, and touch upon some of its interesting features. 
One example feature, and thus potential signature of mod- 
ified gravity, concerns the growth of structure. In Section [3] 
we look deeper at this growth characteristic and attempts 
to generally parameterise it analogous to the equation of 
state for dark energy. We highlight how this extra richness 
in modified gravity can break the observational degeneracy 
with dark energy models and discuss the ensuing limita- 
tions. In addition, we note how modified gravity alters the 
relationship, relative to GR, between the power spectrum 
of the potentials and the matter power spectrum which is 
implicit within weak lensing. Section|4]introduces weak grav- 
itational lensing and its particular importance to modified 
gravity given that it is sensitive to the expansion history, the 
growth of structure and the power spectrum. In this section 
we al so detail t he su rvey and data that we use (CFHTLS- 
wide: IFu et al.l (|2008h - From here on F08) and follow with a 
discussion of the working caveats, including non-linearities, 
and how this data in particular circumvents the issue. Sec- 
tion [5Tl] subsequently contains the analysis and constraints 
of the general braneworld-like model and parameterisation 
of growth through lensing. It is promptly followed in Sec- 
tion 15.21 by the addition of BAO and Supernovae data to 
improve upon these constraints and break parameter degen- 
eracies. For all the analyses in Section 15.11 and Section 15.21 
we implement a Monte C arlo Markov Chain (M CMC) ap- 
proach with CosmoMcH (|Lewis fc Bridle! 120021 ') where the 
resulting plots have been produced with CosmoloGUfl. In 
Section [5.31 we briefly highlight potential systematics in the 
data and quantify any effect on the constraints. We also look 
beyond present day constraints on gravity in Section [6] and 
see h ow the highly exciting futu re weak lensing prob e Eu- 
clid (|Refregier et al.l (|2008l ) and ICimatti et al.l (|2008l )') will 
be able to distinguish between General Relativity and other 
models of cosmic acceleration and growth. We finish in Sec- 
tion[7]with a summary of the paper including a discussion of 
the caveats and limitations as well as suggestions for future 
work. 



^ http:/ /cosmologist.info/cosmomc/ 

^ http:/ /www. sarahbridle.net/cosmologui/ 



2 MODIFIED GRAVITY 

General relativity itself is a modification of gravity. It su- 
perseded the previous established theory of gravity, New- 
ton's Law of Gravitation, with a breathtaking physical prin- 
ciple and mechanism for gravitation phenomena. Although 
an enormous and elegant change in how we think about 
gravity it was quite simply necessary: The previous gravita- 
tional framework did not explain all gravitational processes. 
For example, it did not account for the anomalous preces- 
sion of Mercury. Given the success of Newton's theory many 
attempts were made to understand this within the existing 
framework. In fact, even a form of dark matter was invoked 
in the form of an unseen planet to cause the required pro- 
cession. 

Today one may argue we face a similar choice with the 
evidence of accelerated expansion in the Universe. Again 
early attempts have tried to incorporate some new mat- 
ter component within the formalism of our current the- 
ory of gravity. The simplest procedure has been to intro- 
duce a cosmological constant -perhaps arising from vacuum 
energy-to the usual Einstein field equations. However, a 
disagreement of 120 orders of magnitude represents a se- 
vere fine tuning problem. Other similar avenues have in- 
cluded the introduction of a dynamical scalar field which 
is either trapped withi n a false vacuum or slowly rolling 



down a potential (E.g . IWetterichl (llQSSi'). [Peebles fc Ratr j 
(198^), 'F rieman et "aiT (119951 ). iFerreira fc Jovcel (| 19981 ) and 
[Albrccht fc SkordisI (|2000l )). These Quintessence or dark en- 
ergy models can potentially lead to the desired acceleration. 

Alternatively the more radical, but historically success- 
ful, route is with another modification to gravity. Starting 
from the assumption that any viable theory is described by a 
Langrangian one might at first consider adding terms to the 
Ricci scalar (R) in the Einstein-Hilbert action for General 
Relativity given below. 
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In fact this procedur e was performed so me time ago in 
the context of inflation (jStarobinskiil Il98[]h but has more 
recently bee n analysed for the la te-time low curvature uni- 
verse, in e.g. I Carroll et al.l (|2004l ). where the term 1/R was 
added. It was found to have the desired effect of acceler- 
ation but is ultimately unfeasible as a realistic alternative 
due to its fa i lure t o comply with solar system constraints 
(jChiba et all l200d ). A generalised modification could as- 
sume that the correct underlying theory must be some gen- 
eral function of the Ricci scalar. These models, called f(R) 
models, are being studied extensively in the literature. One 
could potentially generalise this even further to functions of 
the Ricci tensor i?^^ and curvature tensor R^ivpa, resulting 
in f(R,7iji^i?'"',_R^i,p<ji?'"''"') gravity. However, this general 
gravity suffers from h igher order inst abilities, through Os- 
trogradski's theorem (|Woodardll2007l ). and so analysis has 
tended to focus on f(R). It is not exclusively this subset 
of theoretical space that suffers theoretical problems how- 
ever. Healthly theories of gravity seem to be particularly 
rare with most suffering from a whole host of theoretical 
afflictions, from ghost negative ener gy states to tachyonic 
behaviour (|Durrer fc Maartensll2007l ). It is also worth not- 
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ing that if any of these remaining problem-free theories are 
successful in explaining the observed acceleration we are also 
faced with the challenge of explaining why vacuum energy, 
mentioned earlier, does not gravitate! 

Beyond generalising the function in the Lagrangian we 
could also look to higher dimensional models. Within the 
context of cosmology this braneworld scenario can somewhat 
be described as string theory inspired. For example, normal 
matter might be confined to a 4-dimensional brane, where 
the conservation equation p + 3H{p + p) = holds firm, 
but gravity is free to roam into a higher dimensional bulk. 
One pa rticular braneworld scenar io, the Randall-Sundrum 
model (|Randall fc Sundrum|[l999h . approximates to stan- 
dard gravity at low energies but becomes five dimensional 
for high energy and small length scales. For late time accel- 
eration however we desire something that will change over 
large distances and low en ergy scales. The P GP modelf] of 
Dvali-Gabadadze-Porrati ([Pvali et al.l l200oh . described by 
the Lagrangian in Equation is exactly this. 



Ig = — ^ [- / d^-y^J?(^) + f A\^R 

iDTTLr ITc Jbulk Jbrane 



(2) 

It was originally created consisting of a 4-dimensional 
Minkowski brane within a 5-dimensional Minkowski bulk 
and with no motivation t owards explain ing dark energy. 
However the generalisation (|Deffavetll200ll ) to a Friedmann- 
Robertson-Walker brane gave rise to a self-accelerating solu- 
tion. Gravity leaking from this 4D brane into the bulk over 
large scales gives rise to the acceleration through a weaken- 
ing effect. The resulting Friedmann equation is a modifica- 
tion to the General Relativistic case and is given by Equation 
lISj with Tc, the cross over scale, specified in Equation Q. 



H 



2 H SttGp 



(3) 



(4) 



With this modification we have a full description of the 
expansion history. This also allows us to work towards un- 
derstanding the growth of large scale structure giving two 
observational signatures that would enable us to carry out a 
cosmological study. The difference in the background accel- 
eration itself is enough to produce a difference in the growth 
of structure. This can be seen in the second term, the Hub- 
ble drag, in the growth of density perturbations S, for GR 
below, 



5 + 2HS = 47rGp™5. 



(5) 



However, assuming that the only alteration is via 
changes in H is fortunately incorrect. It is fortunate be- 
cause it is the extra modification that allows us to break the 
degeneracy between some general dark energy within GR, 
which can replicate any desired expansion history, and this 
modified gravity model (this will be addressed in more detail 
in Section |3]). The correct approach regarding the evolution 

^ See iLud l l200(i ) for an extensive review. 



of perturbations in this gravitati onal framework is particu- 
lar lY_diffi£uIt_aBd_w^ (|2004h and then 
bv lKovama &c Maartend (|2006l ). It was found that treating 
gravity as 4-dimensional, which leads to Equation ([5]), in- 
duces an inconsistency in the 4-dimensional Bianchi iden- 
tities. Instead with the full and complete five-dimensional 
analysis, and assumptions of a quasi-static regime and sub- 
horizon scales, they found the metric perturbations on the 
brane to be. 
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It is this P factor within Equation ([7} that breaks the 
expansion degeneracy and modifies Equation ([5)l becoming 
rather, 



6 + 2HS = 4nG{ I + 



3p)' 



(9) 



One can see the effect of this analysis by looking at Fig- 
ure [1] We have plotted the linear growth factor for LCDM, 
DGP and a dark energy model with the same expansion 
history as DGP. It is evident that the expansion history has 
considerable infiuence on the linear growth of structure with 
a suppression in the dark energy model relative to a cosmo- 
logical constant. The efi'ect of the five-dimensional pertur- 
bation analysis adds to this suppression and acts to clarify 
the deviation between the dark energy and DGP model. 

With the modified Friedman equation, the correct linear 
growth equation and several caveats it is now possible for one 
to perform tests on the expansion history and/or large scale 
structure for this particular modification to gravity. Some 
of these tests already exist and it has been found that DGP 
is un der tension f r om th e recent influx of cosmological data 
(E.g. ISong et al.l |20o3) and references therein). It is also 
worth noting that this model is potentially not without some 
of the th eoretical proble ms a lluded to above with no tions of 
a ghost jKovamal (120051 ) andlGorbunov et al.l l|2006l )') and a 
strong coupling problem (|Rubakovl 20031 ). 

We do, however, go beyond DGP as an isolated theory 
and examine a phenomenological model that is motivated by 
the concept of an extra di mension with infi n ite ex tent. This 
model, first introduced bv iDvali fc Turner] (|2003l ). interpo- 
lates between LCDM and DGP and allows for the presence 
of the extra dimension through corrections to the Friedmann 
equation with the addition of the parameter a shown in 
Equation (|10|l and rc in Equation (lllf) . 
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where = (1 - fl„,)^ Ho\ (11) 

It is clear that in this case LCDM is recovered when a = 
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Figure 1. g(a) = 5(a)/a, the linear growth, is plotted for a range 
of late time acceleration models. The solid line demonstrates the 
growth for LCDM, the dashed for the five-dimensional braneworld 
model DGP and the dotted for a dark energy model with iden- 
tical expansion history to DGP (uiq = —0.78 and Wa = 0.32 
where w{a) = wq -|- (1 — a)wa)- The difference in the expansion 
history seems to give a significant suppression in growth rela- 
tive to a pure cosmological constant. The effect, however, of the 
five-dimensional perturbation analysis not only adds to the sup- 
pression for DGP but breaks the degeneracy between itself and 
the attempt at replication from the smooth dark energy model. 
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Figure 2. g(a) = S{a)/a, the linear growth, is plotted for var- 
ious values of a that characterise the phenomenological general 
braneworld-likc (mDGP) model. The solid line demonstrates the 
growth for a = (LCDM), the dashed for a = 0.25, the dash- 
dotted for a = 0.5 and the dotted for a = 1 (DGP). Once again 
it is evident that the more DGP-like end of the a spectrum expe- 
riences more suppression in the growth of density perturbations. 



and DGP when a = 1. Furthermore, it is worth noting 
that a < leads to effective equation of states less than 
— 1, whereas a > 1 acts to disrupt both the long matter 
era needed for structure formation and the limits set by 
Big Bang Nucleosynthesis (BBN), while a > 2 corresponds 
to early universe braneworld modifications. Using Equation 
(|10|) it is possible to test the model with probes of the expan- 
sion history and i ndeed this has already been performed by 
lYamamoto et all (|2006l ) with Baryon Acoustic Oscillations. 
If one wants to go beyond this and include tests of large 
scale structure then one needs a formalism for the growth of 
density perturbations analogous to Equations (0, (jS)) and 
therefore ([9]). The problem in this scenario is that in order 
to deduce the growth of perturbations one needs an under- 
lying covariant theory and all that exists i n this modified 
DGP model (mDGP) is a parameterisation. iKovamj (|2006l ) 
circumvented this obstacle by taking DGP (a limit in the 
model) as the structure of the theory. It was subsequently 
found that the metric perturbations take the same form as 
Equations ((6)1 and ^ but instead with, 

Figure [2] demonstrates how the growth of density pertu- 
bations alter within the mDGP model - from LCDM (a = 0) 
to DGP {a — 1). As in the previous Figure it is clear that 



there is a suppression of growth at the more DGP end of 
the a spectrum. 

Although a parameterisation of a general large extra di- 
mensional model the mDGP model now has a definite Fried- 
mann equation that governs the expansion history and, as- 
suming the structure of any underlying theory, it has a set 
of metric perturbation equations and a corresponding den- 
sity perturbation equation. We can therefore treat this as 
a specific model which we choose to constrain later in the 
paper. It is worth noting that using this as a measure of 
deviation from General Relativity or as a parameterisation 
of general modified gravity is not our aim. This would con- 
stitute a poor choice of parameter given the severe lack of 
generalness and incompleteness with regards to the concept 
of modified gravity as a whole. We do however touch upon 
the idea of parameterising modified gravity in the next sec- 
tion. This model has, however, been extremely illustrative 
with regards to the extra richness that can occur in modi- 
fied gravity. Not only does it have varying expansion histo- 
ries but potentially a whole range of perturbation equations 
which alters the growth of structure and the relationship 
to the power spectrum. This is particularly useful when at- 
tempting to distinguish between LCDM, general dark energy 
and modified gravity, and insightful to the probes that will 
be most adept at detecting them. Again, we discuss these 
issues in the following sections. 
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3 GROWTH 



The alteration in the growth of structure within the mDGP 
model demonstrated an additional observational character- 
istic that allows us to further constrain the model and po- 
tentially break the degeneracy with an expansion replicating 
dark energy. It also highlights the possibility of searching for 
signatures of modified gravity in current data by looking for 
changes in the growth of structure. However, even if we do 
have a plethora of new and viable physical principles that 
explain dark energy naturally we might want a way to quan- 
tify how these models affect this growth signature. If we do 
not have the models we might just want to test the influx of 
data for hints of something unusual. Either way it is possible 
to parameterise this extra growth. 

This notion of parameterising growth is analogous to 
the familiar parameterisation of the background expansion 
into Wo and Wa- This is sufficient in describing and restrict- 
ing the multitude of possible dark energy models and ex- 
pansion histories. It is now common procedure to examine 
data and convert it into constraints on various cosmolog- 
ical parameters including wq and Wa- We might therefore 
like to extend this parameter space and allow for the ad- 
ditional signatures of gravity. One possible parameterisa- 
tion for growt h is given by 7 in Eq uation (|13p and was first 
introduced by iPeeblesI (Il980l) and iLahav et alj I (Il99lh and 
later discussed in Wang fc Steinhardtl dlQQSh. Lindej lj2005[ ). 
iHuterer fc LindeJ ifioOTh and lLind"er fc CahnI (|2007l l. 



;5(a)=exp(^^ (^,,(0)^ - 1) dlna) (13) 



By once again looking at Figure [T] we can see that the 
growth factor g(a) is affected by the expansion history and 
by the gravitational framework. It is worth noting that the 7 
parameterisation distinguishes the two contributions to the 
growth, encapsulating the latter in isolation. This is due to 
the impact from the expansion being absorbed into Qrn.(a) 
thus leaving 7 to pick out any remaining remaining con- 
tribution. It is in this way that 7 has become known as a 
modified gravity or beyond-Einstein parameter. It is easy 
to see why given that it detects changes to the growth not 
associated with expansion. This could be down to a change 
in the force law acting on matter represented, for example, 
by the extra factor in Equation Q. And as we alluded to 
earlier, evident in Figure [T] this allows us to distinguish be- 
tween dar k energy and mo d ified gravity. However, as high- 
lighted in iKunz fc Sapond (j2007l l there exists an interest- 
ing caveat. They found that contrary to Figure [l] one could 
force some arbitrary and generic dark energy to replicate 
the growth of DGP. This was achieved by allowing for dark 
energy models with lower sound speeds (c^ 7^ 1) which in 
turn induce clustering in the fiuid. The clustering instigates 
a deepening of the gravitational potential wells thus leading 
to a magnification in the metric perturbations and subse- 
quently an increase in the growth. In addition the existence 
of anisotropic stress was permitted which had the effect of 
suppressing growth. With a careful balance between stress 
and sound speed they succeeded in replicating g(a) for DGP. 
Now although highly fine tuned it is worth keeping in mind 
that observationally detecting some non-LCDM growth fac- 
tor, or 7, would therefore not necessarily constitute modi- 
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Figure 3. g(a) = S{a)/a, the linear growth, is plotted for various 
values of 7 the parameterisation of growth resulting potentially 
from a change in force law. The solid line represents the growth 
for the LCDM paradigm [wq = —1, Wa = 0) with the corre- 
sponding growth parameter of 7 = 0.55. The dashed line shows 
the growth for 7 = 0.68 which is the same as resulting from fiat 
DGP but with the same expansion as in LCDM. The dotted line 
also has 7 = 0.68 but now with wq = —0.78 and Wa = 0.32 thus 
completely specifying the growth of the example DGP model. Fi- 
nally the dot-dashed line shows the growth for a LCDM expansion 
but 7 = 0.45. It is quite clear that a high value of the growth pa- 
rameter corresponds to a suppression of growth. This potentially 
arises from a weakening of gravity. 



fied gravity. In this way, unless one allows for only smooth 
non-clustering dark energy the growth parameter is not just 
a modified gravity parameter. Instead it has the ability to 
pick up on clustered dark energy and modified gravity both 
of which would be interesting. Given the above it is therefore 
our intention to test the growth parameter and see whether 
current data or a future probe can pick out this subtle but 
potentially important eff'ect. Indeed there exists a few early 
attempts, including constraints from peculi ar velocity mea- 



surem ents from low redshift Supernovae |Abatc & Lahav 
l2008f) a s well as future survey for ecast s from|Aincndo lg. ct al.l 
(|2007l ). lHuterer fc Lindeil (|2007D and lHeavens et all (j2003). 

Figure [3] demonstrates the result of varying this param- 
eter on the linear growth factor. The growth for standard 
LCDM corresponds to 7 = 0.55 whereas for flat DGP it cor- 
responds to 7 = 0.68. In this way it is clear that a higher 
growth parameter results in a suppression of growth. 

It is worthwhile noting that other attempts at pa- 
rameterising modified gravity have been made which 
aspire to encapsulate the properties of gravity simi- 
lar to the Parameterised Post-Newtonian (PPN ) param- 
eters for local gravity constraints (|Willl Il993l ). For ex- 
ample, these include parameterising the relationship be- 
tween the two metric potentials ((/> and tp) and/or quan- 
tifying any modification to the Poisson equation (E.g. 
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lAmendola et al. I (l2007l) . IHu fc Sawickil hOOlh . 



II2006I'). I Jain fc Zhang] (l^jUUYI 
iBertschinger fc ZukinI (|2008h l 



( 20071 ). iDaniel et al. 



Ishak et al.l 



( 20081 ) 



and 

In fact these observations 
lead us on to the final modified gravity signature and charac- 
teristic we must consider before our analysis with weak lens- 
ing. For any deviation in the Poisson equation or between 
the metric potentials causes a deviation in the relationship 
between the power spectrum of the potentials and the power 
spectrum of the density contrast. Failure to account for this 
by not modifying the corresponding lensing equations will 
render any ana lysis incomplete. To un derstand this we use 
the notation of lAmendola et all ()2007l ) but in doing so note 
the equivalence of their Q and ri (defined belo w) to GsEf and 
n give n in the other thorough consideration bv lJain fc Zhan j 
i2007l ). 

Firstly, as will be further detailed in the following sec- 
tion, let us notice that in weak lensing, our probe of choice, 
the deflection of light is sensitive to the sum of the metric 
potentials (p + tp. Therefore, we require the power spectrum 
within the lensing statistic (Equation (|23|l ) to be the power 
spectrum of (j> + tp, written P(0_|_^), which we then relate to 
the matter power spectrum Ps. Defining the matter power 
spectrum in Equation (|16|) and the potential power spectrum 
similarly it is obvious that a general relationship between Pg 
and P(^+^) relies on the relationship between 4> + ip a-nd 5 
which in turn depends on the Poisson equation and the re- 
lationship between (p and xp. This is where the parameters 
Q and r] (or equivalently Gcff and vi) are particularly use- 
ful. Here Q parameterises any modification in the Poisson 
equation relating the metric variable (p to the matter den- 
sity p. rj on the other hand, also defined below, describes the 
relationship between (p and 'ip. 



k^fp = -~A-KGQpa5 



%P = {l + ri)4, 



(14) 
(15) 



Now adding the two metric potentials together, substituting 
ip for 77 and (p and transforming the combined Poisson equa- 
tion, as given in the definition of the matter power spectrum 
below 



< 5{k^,z)5{k2, z) >= {2TvfS{ki + k2)Psik, z), (16) 

and similarly for the potential, it is possible to see the gen- 
eral relationship between the power spectra. 
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S7rG)Va"[Q(l + f)]^P,(fc,z) 
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(17) 



W e can, following the notation of lAmendola et al.l 
([2003), define E = Q(l + -1), giving the modification to 
the power spectrum more succinctly as. 



(<t>+i') 



{k,z) 



iTvGfp^a'^T,^Ps{k,z) 
k-i 



(18) 



In a standard cosmological scenario, such as LCDM for 
example, one will have rj = and Q — 1 rendering E = 1. 
This results in the standard relation between the power spec- 
tra as often assumed in the literature. It is clear therefore 
that neglecting E is tantamount to constraining the subset 



of modified gravity models that do not alter the power spec- 
trum relation from the General Relativistic case. Surpris- 
ingly the mDGP model studied earlier is one such model. 
This can be seen by simply adding Equations ([6} and ((Tj 
together and observing the cancellation in /3. Therefore, for 
this model E = 1 also and as a result just the expansion 
history and growth will be altered. Generally however if one 
strives to include general models E should be allowed to 
vary. Note that it acts to modify the amplitude of the power 
spectrum and so a constant value is degenera t e with erg. 
Accordingly, as introduced in lAmendola et al.l (|2007l ). we 
consider the general parameterisation given by. 



E(a) = 1 -I- Eoc 



(19) 



It is our intention therefore to constrain and forecast 
for the specific mDGP gravity model and then, separately, 
constrain the general characteristics of modified gravity. For 
the latter we choose to set E = 1 (or Eq = 0) due to limi- 
tations in data and constrain wq and 7 signatures only. We 
later include Eo 7^ for the Euclid forecasts. It is worthwhile 
noting that 7 and E are functions of the more fundamental 
Q and rj and that these par ameters themselves are , most 
generally, Q{k,a) and ri{k,a) (|Amendola et al.]|2007h . 



4 WEAK LENSING AS A COSMOLOGICAL 
PROBE 

The deflection of light by mass is given by the transverse 
gradient of the metric potentials integrated along the path 
length. 



J d{^ + <P)ds. 



(20) 



This acts to not only change the apparent position of some 
point source but, in turn, distort the shape of distant source 
galaxies. Indeed, one can relate the observed position of the 
image 9i to the true position of the source Os, in the plane 
of the sky, by 



9i = es + 



DjXs ~ X ) 
D{X) 



(21) 



where D{x) is given by the comoving angular diameter dis- 
tance. The subsequent image distortion is given by the dif- 
ferential of this lens equation resulting in the Jacobian 



A 



1 — K — 71 
72 



72 

1 — K -I- 71 



where the convergence n is reconstructed from the shear 7 
(where 7 = 71 -I - ^72 ) which in t urn is measured from galac- 
tic ellipticities (|Refregieil l2003t ) . From the point of view of 
a cosmological model it is the convergence that we are in- 
terested in. This convergence, which represents a weighted 
pro jected mass distribution on the plane of the sky, is given 
(B artelmann fc Schneideill200ll ) for a general mass distribu- 
tion by 



2c2 



Xs 



X) 



(l + z)5(x). (22) 
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Given that we wish to analyse this distribution in a sta- 
tistical way we use the defintion of the power spectrum - 
analogous to Equation (|16p - to find the expression for the 
convergence power spectrum 



4c4 



La 



(x)J 



Ps 



X 



(23) 



where the geometric comoving angular diameter distance 
terms have been absorbed into g{x)- It is therefore now clear, 
given our earlier discussion, how weak lensing as a cosmolog- 
ical probe is particularly useful in studies of modified gravity. 
For this statistic is sensitive to the growth of structure via 
the presence of the linear growth factor squared (s^(a)) in 
the matter power spectrum Pg. It is sensitive to the expan- 
sion history through the terms in the square brackets and 
through the Hubble drag in the growth of structure and also, 
as discussed at the end of the last section, it is sensitive to 
the relation between the power spectrum of the potentials 
and density. In the equation above the relation from Ptf,+^ to 
Pg has already been performed assuming GR as given rou- 
tinely in the literature. Again, it is worth reiterating that 
if there is a modification to the Poisson equation and/or 
to the anisotropic stress one must augment this powerspec- 
trum with the approapriate prefactor given, for example, as 
in Equation (|17p . In addition to these sensitivities it is worth 
adding that because the defiection is given by the potentials, 
which are sourced by mass irrespective of being baryonic or 
dark, weak lensing does not suffer from any unknown bias. 
That is, it probes the entirety of the mass distribution. 

While this probe, in principle, is excellent for our 
chosen study it is worthwhile noting that the shear sig- 
nal is a small 1% distortion on the already existing in- 
trinsic ellipticity. This provides a thorough technical chal- 
lenge that is being combated with a combination of large 
galaxy num ber analyses an d refin e d shear measu r ement 
tech niques (iHevmans et al.l (|2006l ). iMassev et"al] (120071 ) 



and iBridle et al] (|2008l ')'). Further stiU, the first detec- 
tions of weak lensing are p ar ticularly recent (|Bacon et al.l 
(|2000D. iKaiser et al.l (|2000D . IWittman et"al] (|2000l ) and 



Ivan Waerbeke et al. (|200Cl )') and so in this way lensing is 
very much a highly promising, yet developing, cosmological 
probe. Despite this there are already a number of papers 
in the literature that have addressed the relationship be- 
tween weak lensing and mo d ified gr avity or dark energ- 
such as lUzan fc Bernardeaul (l200lf ). ISchimd et al.l i200. 



iDore et al.l (120071'). ISchimd et aLltoOTl'), lAmendola et al 
(120071' ) ■ iJainfc Zhand (|2007t ) and iTsuiikawa fc Tateka~ 
20081 ). With these studies and potential modified gravity 



attributes in lensing it is imperitive to realise that there 
does exist a severe caveat. This is due to the fact that lens- 
ing probes into the non-linear regim e. We are f ortunate to 
be able to use a fitting function (E.g. ISmith e~a l. 2003) for 
the non-linearities in standard gravity. Unfortunately this is 
poorly understood in any deviation from the current frame- 
work and subsequent implementation of the fit would, tech- 
nically, be invalid. Therefore until N-body simulations have 
been undertaken that could generalise the fitting function or 
accurately quantify deviations from it we must strive to work 
in the linear regime as much as possible. Although some at- 
tempts have been made at quan ti fying the v alidity of the 
present fits (E.g. iLaszlo fc BeanI (|2008l ) and lOvaizu et al.l 



(|2008l )) we keep, for now, a strict and linear only analysis. 
There are, in addition, other benefits in avoiding the inclu- 
sion of small scales such as the presence of intrinsic elliptic- 
ity correlations, shear-shape correlations and the presence of 
non-Gaussianity in the error. We therefore utilise the data 
provided by F08 based on the CFHTLS-wide survey which, 
due to its range of large angular scales (up to 230 arcmin- 
utes) probing the more linear regime, is ideal for work on 
non-LCDM cosmology such as this. 



4.1 CFHTLS 

The Canada-France-Hawaii Telescope Legacy Survey Q 
(CFHTLS), based on the MEGAPRIME/MEGACAM in- 
strument, is an ongoing survey with a target of 450 nights ex- 
tendin g over 5 years. The recent analysis bv lBeniamin et ali 
(|2007h has gone beyond t he in itial releases and investiga - 
tions bv lSemboloni et all (|2006l ) and lHoekstra et al.1 (|2006l ) 
which themselves were successful in deriving constraints 
on the ilm — as degeneracy and demonstrating the evolu- 
ti on of the shear signal with redshift. This was achieved 
m iBeni amin et"al] (|2007l ) through a better understanding of 
the redshift distribution and having an increased area. This, 
while marking significant progress, is still not the most op- 
timal lensing analysis for this work. This is because they 
are potentially sensitive to the growth of structures on non- 
linear scales which, as we emphasised above, is undesirable 
for a current study of beyond-Einstein cosmology and weak 
lensing. We therefore loo k to the 3rd ye ar CFHTLS-wide 
release (T0003) given by IFu et al.| (|2008|) (FOB). Ahhou gh 
having a smaller field of view than Benjamin et al.l (|2007l ) it 
utilises much larger angular scales (into the linear regime) 
also avoiding many of the potential systematics mentioned 
at the end of the last section. It is because of this that 
both works reveal approximately equivalent cosmological 
constraints and little constraining power is lost. 

The current sky coverage of 57deg^, approximately 35% 
of the final CFHTLS target area, is reduced to 34.2deg^ af- 
ter masking and the removal of various contaminants. Even- 
tually including five bands this i' band study stretches to 
a magnitude of i'j^g = 24.5 and encapsulating nearly 1.7 
million galaxies has an effective galaxy number density of 
n = 13.3gal/arcmin^ . The data (F08) comes in the form 
of several two point statistics which are relevant to this 
study. We choose to utilise the E correlation function which 
is shown in Equation (|24p and displayed along with the cos- 
mological best fit in Figure |4] 



iE = — I lP4l)MW)dl 



(24) 



As for the aperture mass < Mfp > and shear top hat 
variance < |7|^ > two point statistics this is a weighted 
transform of the convergence power spectrum. In this case 
it is given by a zeroth order Bessel function of the first kind 
Jo . It is in this way that the two point functions vary in their 
sensitivity to various aspects of the power spectrum and in 
turn any systematics. suffers from a constant offset re- 
sulting from a mixing of E and B-modes. A finite survey 



http: / www.cfht . hawaii . edu / Science/CFHTLS / 
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size introduces a maximum angular scale which prevents a 
complete calculation of the shear correlation function over 
larger ranges. This is needed for a separation of E and B 
l|Kilbinger et al.ll2006l ). To alleviate this we alter the statis- 
tic to + c' including the constant offset c' as an ex- 
tra parameter. An expression is then obtained for the offset 
which represents the best fit offset (dx^/dc' = 0) for each 
parameter c hoice. This constitute s an analytic marginalisa- 
tion over c' (|Lewis fc BridldbOOd 'l. We subsequently find the 
expression to be 



(25) 



where the element is the model correlation function, Di 
the data and C the full covariance matrix between all the 
elements. Furthermore, the B correlation function which 
describes the curl component of the shear field, as opposed to 
which measures the curl-free component, is expected to 
be non-zero only for no n-lensing contributions to the shear 
l|Crittenden et al]|2002l '). ft is because of this that is an 
excellent check on any contamination of the lensing signal. 
F08 found no real B-mode contribution except for the pres- 
ence of a very small signal at large angular scales. They find 
however that their cosmological conclusions are not affected 
by this potential mode. In addition, it was shown in F08 
that there was no significant deviation in cosmological con- 
strai nts across any of the aforementioned two point statis- 
tics. iDore et al.1 (|2007l l. who also looked at modified gravity 
models in the context of this CFHTLS-wide data, came to a 
similar conclusion. It is worth noting that other cosmological 
studies of this da ta set includ e a phcnomenological modified 
gravity analysis (|Daniel et al. 200 8) and more recent ly an 
early study of the neutrino mass (fTereno et al.ll200j ). We 
therefore choose for simply to use only one statistic (^b). 

The redshift distribution of the source galaxies, which 
weak lensing is critically se nsitive to, has bee n calibrated 
for the CFHTLS study using lllbert et al.l (120061 ). We decide 
to follow F08 and model this redshift distribution using the 
function 



n{z) = A- 



where A : 



z'' + c 



■dz 



(26) 

where A is the normalisation and a, b and c are three ex- 
tra parameters to be varied and marginalised over in the 
cosmological fit. They found that Equation (|26|l enables a 
closer fit to the distribution data that other common n(z) 
fitting formulae. The observed normalised redshift distribu- 
tion, along with the fitting function evaluated at the best fit 
points found in the mDGP combined probes cosmological 
run (Section 15. 2|) . is shown in Figure [S] 



5 CONSTRAINTS 
5.1 Lensing 

Having earlier discussed the characteristics of certain late- 
time acceleration models in Section |2] and Section [S] and 
having chosen a probe that is potentially capable of picking 



xlO 



Data 
--■BestFit 







9 (arcminutes) 

Figure 4. The open circles with associated errorbars represent 
the ^E two-point statistic as a function of 6 (arcminutes) for the 
CFHTLS-wide survey used in this paper. We however selectively 
use scales greater than 30 arcminutes to remove the unknown 
non-linear effects and in doing so strive towards a more linear 
analysis. The red dashed line shows the best fit values as found 
with the combined probes mDGP analysis (Section 15.21 1. 



out these particular behaviours with weak lensing, we are 
now in a position to make cosmological constraints based on 
the data which we decided, in Section [4. II was the most suit- 
able. That is, we start by making an analysis of the mDGP 
gravitational model and then separately the parameterisa- 
tion of growth with the F08 CFHTLS-wide lensing data. 

We perform a full likelihood analysis using a Monte 
Carlo Markov Chain (MCMC) approach on a set of 7 cos- 
mological parameters for the mDGP analysis and 8 for the 
growth parameterisation analysis. We vary Qrn, h, erg, a, 
b, c which are common to both models, in addition to a 
for mDGP and wq and 7 for the growth. The wo is not 
included for the former model as any given a uniquely spec- 
ifies its own expansion history as it does for the evolution of 
perturbations-from GR to DGP. We decide not to vary the 
spectral index Us because the data needed at large angular 
separations will be insufficient for any constraint. Instead 
we set it to Us — 0.963 consistent wit h the best fit five-year 
WMAP result (jDunklev et al.1 12OO8I ). Likewise we neglect 
varying Wa and Eq for the growth model due to limitations 
in the data. Instead we set them to Wa ~ and Eo = 0. 
Finally, we assume a flat universe throughout. 

For the lensing analysis the Gaussian log-likelihood is 
given by 



(27) 



where the data vector D is given by the measured ^E{9i)- 
The theoretical prediction deduced at the corresponding an- 
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Figure 5. The black dashed line represents the source redshift 
distribution with associated errors in the bins. The red solid line 
is given by the fitting function in Equation I I26I I. The fit is drawn 
for the function evaluated at the best-fit points as deduced by 
the combined probes analysis for mDGP (Section [53J ■ This cor- 
responded to best fit values: a = 0.614 ± 0.034, 6 = 8.11 ± 0.681, 
c = 0.627 ± 0.0610 and A = 0.6462 consistent with F08. 



gular scale 6 is represented by entries in f and, finally, 
is the inverse covariance matrix for the chosen observable. 
At present variation in either a, 6 or c in the redshift dis- 
tribution is detected implicitely through the alteration of 
the model matter power spectrum. However, in implement- 
ing the MCMC approach regions of parameter space could 
potentially be sampled that would correspond to configu- 
rations of a, h and c incompatible with just the knowledge 
of the redshift distribution alone. We therefore follow the 
procedure in F08 and multiply the likelihood above by the 
likelihood of the redshift distribution given by 



2 

Xn{z) 



{rir - n{zi)Y 



(28) 



where n^, the observed number of galaxies in a bin, are 
shown in Figure O n{zi) represent the values of the fit- 
ting function in Equation (|26|l evaluated at the bin cen- 
tred redshifts. While ignoring cross-correlations in the bins 
we include ai which is the error in n;. This error includes 
Poisson noise, sample variance and the associated redshift 
uncertainty. In additio n to the redshift distr ibution we also 
include an HST prior (jPreedman et al.ll200ll 'l on h as given 
by, 



Xhst — 2" 



1 (/i- 0.72)2 



0.082 



(29) 



For the lensing only analysis we cut the data such that 
we use angular scales greater than 30 arcminutes, which we 



Figure 6. The diagram above demonstrates the attempted con- 
straint on a parameterised gravitational model that is motivated 
by the concept of a large extra dimension (mDGP). The contours 
are for Q,m, the matter density, and a the modified gravity pa- 
rameter where 5 other cosmological parameters (/i, o-g, a, b and c) 
have been marginalised. Here only angular scales greater than 30 
arcminutes have been used in order to avoid the non-linear regime. 
The selection of data is from the CFHTLS-widc (F08) survey with 
analysis using the E correlation two point statistic I^e + c'. Note 
that for this mDGP model a = corresponds to LCDM whereas 
o = 1 is equivalent to the DGP braneworld model. 



reiterate, is to avoid the unknown non-linear contribution 
to the lensing constraint. The analysis for the mDGP model 
is shown in Figure [B] It shows the marginalised ilm, and 
a contours where, as shown in Section [S] a parameterises 
a gravitational model with a large extra dimension. The 
mDGP model interpolates between LCDM (a = 0) and the 
DGP braneworld model (a = 1). It is clear therefore that a 
lensing only analysis is presently not capable of constrain- 
ing mDGP-at least in the context of physically more viable 
models (a < 1). Indeed we find a similar difficulty in 7 with 
no constraint on any reasonable physical values in a lensing 
only analysis. Firstly, we should perhaps not be too sur- 
prised as we have used a relatively new cosmological probe 
and have, in neglecting the non-linear scales, used only a 
third of the data points. Tomographic information is not 
yet currently available for this data either which will vastly 
improve information on the expansion history and hence a. 
We have also allowed significant cosmological freedom with 
the variation of 7 parameters. On the other hand, however, 
this does not mean that lensing, even with a current analy- 
sis, is not useful with regards to our late-time acceleration 
models. In order to see the current effectiveness we now look 
at Baryon Acoustic Oscillations, Supernovae and lensing in 
combination in the next section. Then in Section |S] we see 
how the future weak lensing survey Euclid will improve upon 
today's lensing only constraining power. 
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5.2 Supernovae and Baryon Acoustic Oscillations 

In the previous section we performed a preliminary analysis 
based on linear to quasi-linear weak lensing data alone. This 
probe, while having characteristics significant for the dis- 
crimation of late time acceleration, was unable to constrain 
either the mDGP model or the growth parameterisation. 
As such it is desirable for us to combine it with other cos- 
mological probes in order to improve potential constraints. 
Moreover, it is most beneficial to combine weak lensing, an 
indicator of growth and expansion, with distance indicators. 
This is because of the particular degeneracies that exist be- 
tween the parameters following an isolated study of lensing. 
For example, there exists a degeneracy between w and 7 
and so tighter constraints on the expansion history will act 
to aid any constraint on 7. Furthermore, inclusion of addi- 
tional expansion data will aid the constraints of the mDGP 
model given that different a corresponds to different late 
time accelerations. We therefore choose to include both Su- 
pernovae and Baryon Acoustic Oscillations (BAOs) which, 
due to their vastly different ^Im — w degeneracies, are also 
extremely complementary to one another. 

BAOs, used as standard rulers and observed in the 
galaxy distribution, aim to test cosmology through the 
distance-redsh i ft rela tion. Using the data and notation of 
IPercival et ai] (|2007l ') we look to utilise the distance mea- 
sure given by 

Dv{z) = [{l + zfD\cz/H{z)]^, (30) 

where Da is the angular diameter distance and H{z) is the 
Hubble parameter. Specifically, it is the ratio Ts/ Dv{z) that 
we examine where is the comov ing sound horizon at re- 
combination. IPercival et aD (|2007l l detects the BAO in the 
clustering of 2dFGRS and SDSS galaxy samples and the 
clustering of SDSS LRGs to quantify this measure ai z = 0.2 
and z = 0.35, respectively. For each likelihood evaluation we 
compare this data to rs/Dv{z) calculated with Dv{z) from 
Equation (|30|l and the comoving sound horizon evaluated 
using the formulae in lEiscnstcin fc Ku (1998 ). 

For the inclusion of Supernovae we use the data pro- 
vided from the firs t year Supernova Legacy Survey (SNLS) 
lAstier et al.ll20o"6h . This data set includes 71 type la Super- 
novae also detected at the Canda-France-Hawaii Telescope. 
Here we use the distance modulus /lo, a measure of the lu- 
minosity distance, as the observable. 

Aio = 51oglO(dL(2))-f25. (31) 

This is given in the log-likelihood (Equation (|32p ) also as ^0, 
where is the observed value, aint is given by the intrinsic 
dispersion of the absolute magnitudes and (j{fiB) by peculiar 
velocity and light curve parameter information. 

2_ (/iB-51og,o(dL(e,^))-25)^ 

With this machinery in place we are now able to perform 
additional tests on the mDGP model. We do not constrain 
the growth parameterisation with these probes in isolation 
as in this format they have no growth information. We do. 



however, attempt to constrain the growth with a combined 
analysis at the end of this section. 

By looking at the top left hand panel of Figure [7] 
we can see that it is possible to place a constraint on a 
with a BAO only analysis. It should be noted that the ap- 
parent disfavouring of DGP (a — 1) is not as promising 
as it first appears for we have, replicating the analysis of 
lYamamoto et all (|2006l ) with BAO only, varied just Qrn. and 
a with Qt and h held fixed at 0.044 and 0.66, respectively. 
The top right hand panel of the same Fig ure demonstrates 
tjie^need for caution as we go beyond the lYamamoto et al.l 
f|200(^ ) analysis. By allowing more cosmological freedom (i.e. 
varying fit, and h) and including the Supernovae data the 
two probes are in fact incapable of disfavouring DGP. We 
can now see how, even at present, the weak lensing data is 
useful for the study of this modified gravity. The Qrn — a 
contours for a joint analysis with all three combined probes 
is displayed in the bottom left hand panel. Once again we 
use only angular scales greater than 30 arcminutes to avoid 
the unknown non-linear regime. For the combined analysis 
it is worth noting that we vary a significant total of 8 cos- 
mological parameters: ilm, h, erg, Slj,, a, b, c and a. It is 
now evident that the addition of the weak lensing analysis 
is beneficial with a mild disfavouring of the DGP end of the 
a spectrum. Indeed, we include in the bottom right hand 
panel the ID probability distribution for a, in the process 
demonstrating that a < 0.58 at the 68% confidence level and 
a < 0.91 at the 95% confidence level. This corresponds to 
a disfavouring of DGP at over 2a. Furthermore, we include 
in the left panel of Figure [S] for interest, the same anal- 
ysis but with all angular scales (1-230 arcminutes). There 
is a noticeable but not significant improvement in the con- 
straint on mDGP. The slight increase in disfavouring leads 
to Q < 0.56 and a < 0.86 at the 68% and 95% confidence 
levels, respectively. Again to reiterate it should be noted 
that this analysis, in using all angular scales, ventures into 
the non- linear regime. 

Having had some success with a combination of the 
three cosmological probes it is worth investigating whether 
they can aid the determination of the far more subtle growth 
parameterisation 7. While the BAO and Supernova will not 
add growth information explicitely they may help reduce 
the parameter degeneracies and cut the parameter space 
favourably. We find however, by looking at the right panel 
of Figure [51 that at present and for meaningful values of 
7 we still have insufficient constraining power. This is un- 
derstandable given that for the mDGP model a contained 
growth and expansion information. In this way the Super- 
novae and the BAO actively contributed to the constraint 
while the lensing constrained it through the comoving di- 
ameter distances D{x) in the convergence power spectrum, 
through expansion terms within the growth via the Hub- 
ble drag and finally through the pure growth contribution 
as seen in the P contribution. Constraining 7 on the other 
hand is equivalent to just changes in /3 and is therefore far 
subtler. However, just because we do not have the current 
data to pick out this effect it should not deter us from pur- 
suing this potential signature of modified gravity. In fact, 
future cosmological probes, such as weak lensing, may well 
be able to extract this contribution and tighten constraints 
on beyond-Einstein cosmology. We look to the future in Sec- 
tion [6l 
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Figure 7. The plot in the top left panel shows the constraint on flm and a. Although appearing to disfavour DGP (o = 1) as in the 
analysis by lYamamoto et al.l l|200d ) the remaining parameters fij, and ho have been fixed at 0.044 and 0.66 respectively. We go beyond 
this in the top right panel which contains constraints given on the same parameters but when Supernovae data is added and ili, and ho 
are allowed to vary. We see now that the la contour is beyond the bounds of the plot and so no constraint can be inferred. The benefit 
of the weak lensing data is seen in the bottom left panel where once again we use angular scales greater than 30 arcminutes from the 
CFHTLS-wide (F08) lensing survey. We also vary Om, h, erg, f2j, a, b, c and a whilst keeping Us = 0.963. With this addition it is evident 
that there is a visible improvement in constraint and that DGP is marginally disfavoured. This is exemplified in the bottom right panel 
where we include the ID marginalised probability distribution (solid line). Wc find that the joint analysis gives constraints on mDGP of 
a < 0.58 and a < 0.91 at the 68% and 95% confidence levels, respectively. The dotted line represents the mean likelihood of the samples. 
Finally, the dashed contours in the bottom left hand panel show that our constraints arc roughly insensitive to any systcmatics in the 
data such as an over or underestimation of the shear at high redshift (Section 15.311 . 



5.3 Accounting for systematics 

Despite the absence of any significant B-modes in the 
CFHTLS data there exists a potential underestimation of 
the shear at high redshift (Kilbinger et al. in prep.). In or- 
der to account for this potential effect we use the model 
introduced in Kilbinger et al. and multiply the redshift dis- 
tribution at high redshift {z > 1) by some constant CQ-where 
Co < 1 constitutes underestimation, co > 1 overestimation 
and Co = 1 no alteration in the shear. Following the afore- 
mentioned paper we vary co and marginalise over the pa- 
rameter in an additional cosmological run, thus accounting 
for any such systematic. We subsequently find little change 



in the constraint on mDGP as shown by the dashed contour 
in the bottom left hand panel of Figure [7] Indeed the 68% 
confidence interval is only shifted from a < 0.58 to a < 0.61 
with systematic treatment. Moreover, with the cut in an- 
gular scales at S = 30 arcminutes it has not been possible 
to meaningfully constrain co itself. Further treatment and 
causes of potential systematics in the data are detailed in 
van Waerbeke et al. in prep, and Kilbinger et al. in prep. 
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Figure 8. The left panel is an analysis of the mDGP model with weak lensing, BAO and Supernovae, as before, but in this case with the 
full range of angular scales (1-230 arcminutes) . There is a slight, but not significant, improvement than in the more linear analysis. Here 
we find a < 0.56 and a < 0.86 at the 68% and 95% confidence levels, respectively, ft should be noted that this analysis includes data 
from the unknown non-linear regime however. The right panel demonstrates the current challenge in constraining the gravitational-as 
opposed to the expansion— contribution to the growth of structure. We find that with current data we are unable to put any bounds on 
reasonable values of the 7 parameter. This plot contains an analysis with weak lensing, Supernovae and BAOs. Implicit in this plot is 
the variation of also h, erg, H;,, a, b, c and wq. Only angular scales greater than 30 arcminutes have been used. 



6 FUTURE PROBES - EUCLID 

Within the foreseeable future the age of precision cosmol- 
ogy looks likely to get ever more precise. The field of weak 
lensing is most definitely no exception and it is perhaps set 
to be one of the most promising areas of development. One 
striki ng reason, am ong others, is the planned Eucl id mis- 
sion (|Refregier et al. (2008J and lCimatti et al.l (|2008l )'). This 
space based project aims to have an effective sky coverage of 
20000 square degrees, probing the galaxy distribution with 
a median redshift Zm ~ 0.9 and giving an effective galaxy 
density of 40 gal/arcmin . Hence it is of great interest to 
see how this future survey will probe the nature of gravity, 
whatever that may be. We therefore undertake a Fisher ma- 
trix analysis with the intention of forecasting, with lensing 
only, how Euclid will be able to constrain the mDGP model 
(a) as well as the general parameterisation for gravity (7 
and E). The Fisher matrix formalism for weak lensing is 
given by 

«.-ES;Cov-'g, (33, 

where C is the weak lensing observable shown by 

C^Al) = Ptj+ <lL > (34) 

and P^j is the power spectrum similar to Equation (|23p but 
with indices i and j denoting tomographic bins. Further still, 
rxi is the average galaxy number per steradian in bin i and 
< 7?^t > 2 is the rms intrinsic shear in each component, which 
is equal here to 0.22. dpi denotes the derivative with respect 



to a parameter p that is varied and Gov is the covariance 
matrix given by, 

Cov[C^,{l),C',,il')] = (2; + A/ [CUl)CUn+Ca{l)C-,{l)]- 

(35) 

For our analysis we assume the redshift distribution 
given in Equation (|36|) where zq — 2m/1.412, a — 2 and 
/3 — 1.5. We allow for five redshift bins with divisions such 
as to give approximately equal galaxy number in each bin. 
The photometric error is also accounted for by using the 
parameterisation az — ap{l + z), with (jp taken to be 0.03. 

n{z) oc z"exp(— (z/zo)'') (36) 

For the base cosmology we vary, about their fiducial val- 
ues, the parameters f2m = 0.3, h = 0.7, as = 0.8, flh ~ 0.05 
and Hs = 1.0, which are common to both analyses, with also 
Mo = —1.0, Wa = 0.0, 7 = 0.55 and So = (see Section [Sjl 
for the general parameterisation, and q = for the mDGP 
model. We therefore vary a full 9 parameters for the general 
modifed gravity parameterisation and 6 parameters for the 
specific braneworld model-remembering that a specifies its 
own Wo, Wa, 7 and So = 0. 

The resulting forecasts for this future project can be 
seen clearly in Figure [5] The left panel shows that Euclid 
will be able to put considerable strain on any braneworld- 
like gravity scenario that resembles the mDGP model. The 
solid (la) and dashed {2a) contours are significantly within 
the Q = 1, or DGP, bound. In fact, this demonstrates that 
Euclid will potentially constrain a to within an error of 0.104 
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at the 68% confidence level. This is in constrast to Figure [6] 
where no constraint was possible with a lensing-only study. 
Indeed, note that for this analysis only contributions from 
Z = 10 to Zmax = 500 were considered such that the deeply 
non-linear regime could be neglected. 

The right panel in Figure |9] again illustrates the ex- 
pected constraining power of Euclid but now with regards to 
general modified gravity. For this general parameterisation 
we performed two runs with one corresponding to 1 contri- 
butions from ? = 10 to /max ~ 500 (red contours) and the 
other with contributions from Z = 10 to Zmax ~ 10000 (black 
contours). Here we see that Euclid will be able to extract 
the growth characteristic thus allowing, potentially, a strong 
cosmological test of our gravitational framework. Indeed, it 
will be able to constrain the fiducial 7 — 0.55 (LCDM) to 
within an error of 0.0446 at la with Zmax = 500 which is 
further tightened to 0.038 with Z^ax = 10000. The other 7 
cosmological parameters {h, as, fib, wo, Wa, ris and Eo) have 
been varied and marginalised over. Again this forecast is in 
contrast to Figure [8] where even the combined probes of lens- 
ing, BAOs and Supernovae were incapable of constraining 
the growth. 

Finally, we include in Figure [10] a marginalised con- 
tour for 7 against Eq which further highlights how modified 
gravity or very generic dark energy signatures can be con- 
strained, consistently, with weak lensing. We find that with 
^max = 500 (red contours) this survey could constrain alter- 
ations in the powerspectrum with an error in Eo of 0.25 at 
la. This parameter is more sensitive to the range of scales 
used than 7, however, with /max = 10000 (black contours) 
confining Eo to within 0.069 of the fiducial value Eo — 0. 

We should emphasise th at the spectroscopic element of 
Euclid (|Cimatti et ahlfioOSl l will also be able to co n strain 
gro wth via redsh i ft spa ce distortions (jPeacock et al.l (|2002l ) 
and lGuzzo et all (|2008l )V 



7 DISCUSSION AND CONCLUSION 

To summarise, we have noted that the surprising, but well 
confirmed, late-time acceleration of the universe could be 
the result of a modification to gravity. We then, in Sec- 
tion [2I reviewed the concept of modified gravity detailing 
in the process a model that is motivated and acts to param- 
eterise a large extra dimension (mDGP). Interpolating be- 
tween LCDM (q = 0) and DGP {a — 1) it can be tested as a 
model in its own right and/or used as an example to demon- 
strate the rich set of observational signatures that are likely 
to arise for a modified gravity model in cosmology. These 
signatures include the expansion history, and rather inter- 
estingly, the growth of structure and a modification to the 
relationship between the potential power spectrum and the 
matter power spectrum. With these characteristics in mind 
we then discussed attempts in the literature to parameterise 
modified gravity in this way. This included a growth param- 
eter 7 and a power spectrum parameter E. 

In Section [3] we introduced weak lensing and related 
its attributes to modified gravity given that it is sensitive to 
the expansion history, the growth of structure and the power 
spectrum. We noted a severe caveat in the use of non-linear 
scales and therefore, in Section [4. II described the approapri- 
ate choice of survey (CFHTLS-wide) and data {6 > 30 ar- 
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Figure 10. The above plot shows the marginalised 7— So forecast 
for a weak lensing only analysis with Euclid. These two parame- 
ters, which could represent modified gravity or generic dark en- 
ergy signatures, demonstrate how this future weak lensing probe 
will potentially place firm constrains on any model of late-time 
acceleration. The black contours correspond to imax = 10000, 
demonstrating an error of 0.069(la") on Eq, whereas the red con- 
tours correspond to imax = 500 giving instead an error of 0.25. In 
both cases the inner and outer contours are la and 2a, respec- 
tively. 



cminutes) used in the cosmological analyses. The subsequent 
lensing only constraints were given in Section [5.11 where we 
showed that one could not yet constrain meaningful values 
of Q or 7 with the current data. We then, by adding it to 
Baryon Acoustic Oscillation and Supernovae data, demon- 
strated that weak lensing was highly beneficial in aiding the 
constraint of mDGP in combination. For without the in- 
clusion of the lensing data the expansion-only probes were 
incapable of constraining a = 1. We found however that the 
combined probes disfavoured the DGP model with over a 
95% confidence level where specifically a < 0.58 at la and 
a < 0.91 at 2a. We then showed that a constraint on the 
subtle, yet potentially important, growth signature is be- 
yond the current weak lensing, BAG and Supernovae data. 
We allowed almost total cosmological freedom in all these 
analyses, varying parameters: Sim, h, as, Qb, a, b and c for 
both models, in addition to a for mDGP and wo and 7 for 
the growth model. Furthermore, we used the two point 
statistic while analytically marginalising over the residual 
offset c'. We also demonstrated in Section [5.31 that our re- 
sults are insensitive to a potential over or underestimation 
of the shear at high redshift. 

Finally in Section [6] we looked towards the future space 
based weak lensing survey, Euclid, and discovered that it 
will have significant ability to distinguish between modified 
gravity and LCDM. We included a forecast for the mDGP 
model finding that even for a lensing only analysis Euclid 
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Figure 9. The left panel displays the ability of Euclid's potential constraining power with regards to the mDGP model in a Icnsing only 
analysis. Here the Icr contour (all solid lines) is well within the a = 1, or DGP, line and so it will be easily distinguishable from LCDM 
{a = 0). Indeed this corresponds to an error of 0.104 on a with ^max = 500 (all red contours) in stark contrast to today's constraint. The 
right panel shows the marginalised contours for the general growth parameterisation. Again, it seems that Euclid will provide excellent 
insight into any potential modified gravity signatures. Specifically it is found that it will be possible to constrain 7 with an error of 0.045 
(1(t). This is tightened further to an error of 0.038 when Zmax = 10000 (black contours). The parameters h, erg, Q^, and ris have been 
varied and marginalised over for both models considered here while in addition wq, Wa and So have been marginalised for the growth 
model. 



could restrict a to within 0.104 of the fiducial a = at 
la, even when the deeply non-linear regime has been re- 
moved (Zmax = 500). In addition, we included a complete and 
consistent forecast for generalised modified gravity demon- 
strating that deviations from a fiducial Eo = of 0.25 at 
the 68% confidence level will be possible with Zmax = 500. 
With ^max = 10000 this gets restricted to AEo = 0.069. It 
will also confine 7 to within 0.045 (/max ~ 500) and 0.038 
(^max = 10000) of the fiducial 7 = 0.55 at la, where a full 9 
cosmological parameters were varied. 

In our analyses with data we have, except as an example 
case, actively removed angular scales less than 30 arcmin- 
utes. This was to remove the contribution from the unknown 
non-linear regime from the constraints. This constitutes re- 
moving vast quantities of attainable data. Not only is this 
clearly not maximising the available information but it could 
be that non-linear physics acts to emphasise any difference in 
gravitational theory. Therefore it would naturally be highly 
beneficial to understand the evolution of these density per- 
turbations within this regime. In order to do this there needs 
to be continuing attempts in the community to study N- 
body simulations with varying gravitational frameworks. 

Even though we have detailed the advantages of weak 
lensing in a modified gravity study, and even though Euclid 
in particular will be deeply insightful it seems, however, 
a collected and coordinated assault on our gravitational 
framework would prove even more advantageous. This 
might exist in the form of a combination of pr obes as, 
for example, discussed in I Jain fc Zhand (|2007l ). where 
ideally the four perturbation variables (j), ^, 5 and 6 axe 



independently targeted. This in principle would allow us 
unprecidented experimental scrutiny on the structure of our 
gravitational theory over large scales. However, to conclude, 
it is clear that to actually solve the tantalising problem 
that grips cosmology there needs to be a paradigm shift 
analogous to those that have revealed the great problems of 
the past. 
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